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Abstract. We investigate the phenomenon of multilayer formation via layer-by-layer deposition of alternat- 
ing charged polyelectrolytes. Using mean-field theory, we find that a strong short-range attraction between 
the two types of polymer chains is essential for the formation of multilayers. For strong enough short-range 
attraction, the adsorbed amount per layer increases (after an initial decrease), and finally it stabilizes in 
the form of a polyelectrolyte multilayer that can be repeated hundreds of times. For weak short-range 
attraction between any two adjacent layers, the adsorbed amount (per added layer) decays as the distance 
from the surface increases, until the chains stop adsorbing altogether. The dependence of the threshold 
value of the short-range attraction as function of the polymer charge fraction and salt concentration is 
calculated. 



PACS. 82.35.Gh Polymers on surfaces; adhesion 
tomers, and plastics 



82.35.Rs Polyelectrolytes - 61.41.-|-e Polymers, elas- 



1 Introduction 

The study of polyelectrolyte (PE) chains interacting with 
charged surfaces has generated a great deal of attention in 
recent years. This interest arises, in part, because of the 
numerous biological and industrial applications. The ad- 
sorption and depletion of polyelectrolytes on charged sur- 
faces have been extensively studied using analytical 
Pll^ll^ll^ and numerical [7IIHlinill0llll| solutions of the non- 
linear mean-field equations, scaling considerations |7||S1 

lattice models |18lll9ll2()ll21j and computer simulations |22l 
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,123,24,25]. 

In recent years, formation of polyelectrolyte multilay- 
; ers has been investigated experimentally ■26"27','28"29"30', 

'EI]D2,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47 . Thes(bpposite to the previously adsorbed PE chains^ and washe'd 
multilayers are composed of alternating positively and neg- again ('wash'). The process of dipping and washing: 'A', wash. 



Fig. 1. A schematic illustration of the layer- by-layer PE de- 
position process. The charged surface is first dipped into a so- 
lution of oppositely charged PE chains (marked 'A'). After the 
chains adsorb onto the charged surface, the surface is removed 
and dipped into a clear water solution (marked 'wash'), remov- 
ing any extra non-adsorbed chains. The surface is then dipped 
into another solution of PEs (marked 'B'), carrying a charge 



atively charged PEs and are constructed via a layer-by- 
layer adsorption of polyelectrolyte chains, shown schemat- 
ically in Fig. 1. The first stage is to dip a charged sur- 
face into a solution of oppositely charged polyelectrolytes 
and salt. After the polyelectrolyte chains adsorb on the 
charged surface, the surface is taken out of the solution 
and washed in a clear water solution. The washed sur- 
face and adsorbed layer are then placed in a solution of 
another polyelectrolyte, of an opposite charge to the first 
PE chain (see Fig. 1), and then washed again. This pro- 



'B', wash, 'A', 
dred layers. 



wash,... can then be repeated for several hun- 
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cess can be repeated for several hundred times and results 
in a PE multilayer build-up 32 . More recent studies |30l 
I31| have shown that these multilayers have interesting and 
potentially useful applications both for planar and spheri- 
cal geometries, leading the way to creation of multilayered 
and hollow spherical capsules. 

Theoretical models for multilayer formation have also 
been considered in recent years [HII11II48II49II5()II51II52| . In 
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previous studies [()ll48| . it was suggested that the inversion 
of the surface charge by polyelectrolyte adsorption occurs 
under the following conditions: (i) sufficiently large salt 
concentration; (ii) theta solvent conditions; and, (iii) weak 
short-range (SR) interactions between the surface and the 
PE chains. The charge inversion, together with complex- 
ation of polyanions and polycations have been used to 
model the multilayer formation |48| . 

In a previous publication we reached different 

conclusions by solving numerically the relevant mean- field 
equations for the PE adsorption. We found that surface 
charge inversion will not occur under a broad range of sys- 
tem parameters. More specifically, it was shown that full 
charge inversion occurs only for strong enough SR surface- 
PE interactions. Related results have been reported in a 
separate study [3], where the full charge inversion was 
found to occur only for strongly hydrophobic PE-backbone. 
Namely, PE chains in poor solvent conditions. The strong 
hydrophobicity in the bulk creates an effective attraction 
to the surface, replacing the bare surface-PE short-range 
attraction. 

In the present work, we study the formation of alter- 
nating charged PE multilayers as a function of polyelec- 
trolyte charge, added salt and SR interaction between the 
PE chains. We find that strong SR (non-electrostatic in 
origin) interactions are necessary for the formation of such 
multilayers, and that the adsorbed charge of the alter- 
nating layers is not necessarily equal, or even close to, 
the initial surface charge. Our analysis shows that the ad- 
sorbed amount (per added layer) in the initially adsorbed 
layers always decreases. If the SR interactions between 
the PE chains are too small to attract another PE layer, 
the multilayer formation will stop after a small number of 
layers. However, if the SR attraction between the alter- 
nating PE layers is significant, the adsorbed amount (per 
added layer) starts to increase back, and then saturates 
and forms a stable multilayer stack. We also show how 
the multilayer formation depends on the solution salinity 
and PE charge fraction. 

In the next section. Sec. El we present the mean field 
equations for multilayer formation and the numerical 
method used to solve them. The numerical results for mul- 
tilayer formation follow in Sec. |3 and their discussion in 
Sec.0] We end with conclusions and suggestions for future 
research in Sec. El 



2 The Mean Field Equations 

Consider an aqueous solution in contact with a bulk reser- 
voir of salt ions, and a dilute bulk concentration of long 
polyelectrolyte chains. The solution is in contact with an 
infinite and planar surface. The surface is oppositely charged 
and attracts the PE chains. The mean-field equations for 
this system were formulated in Refs. J,8,9,1Q,, and are 
repeated here. 

^ = K2sinhC-4^?B/0' (1) 



Y0 = -'^^ + /C'^ + -^<^^ (2) 

4>^ is the local monomer concentration, x the distance from 
the charged surface, C, — eip/kBT the renormalized (di- 
mensionless) electrostatic potential, a the monomer size, 
and / the charge fraction of the PE monomers. The Debye- 

Hiickel length = (BTrZsCsait)""^^^ is the screening length 
for electrostatic interactions in presence of salt ions, and 
Ig = e'^/eksT is the Bjerrum length, which is approxi- 
mately tA for water with e = 80 and at room tempera- 
ture. 

Equation |^ is the Poisson-Boltzmann equation, where 
the salt ions obey the Boltzmann distribution and together 
with the monomer charges they act as charge sources for 
the electrostatic potential. Equation is the Edwards 
equation for the monomer order parameter 0, where the 
chains are subject to an external potential composed of 
an electrostatic potential and an excluded volume inter- 
action between the monomers. Note that in Eq. (0) we in- 
cluded both the second virial term modeled by v and the 
third virial one modeled by w^. In most previous studies 
of PE adsorption the third virial term has been omitted, 
because of the dominance of the electrostatic interactions 
and the second virial term. In the case of multilayer for- 
mation, however, the third virial term becomes significant, 
as is explained in the next section. Finally, we note that 
Eqs. ^ and (O are written for vanishingly small bulk con- 
centration of monomers, and under the assumption that 
the ground state dominance approximation holds. 

In order to model the build-up of multilayers, we note 
that the experimental multilayer build-up is done via a 
layer-by-layer adsorption of cationic and anionic PE chains, 
as is illustrated in Fig. 1 and explained in the introduction. 
During the adsorption of each layer, the chains from pre- 
viously adsorbed layers are believed not to dissolve back 
into the solution. Note that the layer-by-layer build-up is 
not a thermodynamically equilibrium process. Any mod- 
eling of this phenomenon should take into account these 
specific stages. 

In our model, electrostatic and short-range (SR) inter- 
actions with the PE chains of the previous layers are taken 
into account. The SR interactions between the anionic and 
cationic PEs may have several origins. The repulsive SR 
interactions include excluded volume interactions, while 
an attractive SR interaction (beside the electrostatic at- 
traction) arises from polyelectrolyte complexation. We do 
not offer in the present work a detailed explanation for 
the complexation origin. We rather assume its existence, 
which yields an effective SR attractive interaction, and 
investigate under which conditions it will lead to multi- 
layer formation. We assume that the electrostatic attrac- 
tion and ion pairing between the chains in the adsorbed 
multilayer and the adsorbing PE chains allow the adsorb- 
ing chains to penetrate the multilayer. This penetration 
slows down the PE chains dissolution back into the solu- 
tion, and creates an effective SR attraction. This interac- 
tion has a non-equilibrium origin, but for dense layers, it 
should last long enough to allow for multilayer formation. 
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The electrostatic interaction between the PE chains 
is taken into account by adding the PE chain charges in 
the Poisson-Boltzmann treatment, Eq. The SR in- 
teractions between the cationic and anionic PE chains are 
taken into account very simply by adding a different inter- 
action parameter x 7^ in Eq. Q). Hereafter, we assume 
that the charged fraction / of monomers for the negatively 
and positively charged PE chains is the same. As the PE 
adsorption is done layer by layer, the equations governing 
the adsorption of the ith layer are: 

d^C 

-j-j = sinh C - 4:TTlBf {ziSs + Zi+iSo) (3) 

~^~dx^ = vSs(t>i + fziC(t>i ~ xSo<l>i + (^'^S'^(j)i (4) 

where C(a;) is the electrostatic potential, and <pf and fz^ 
denote the monomer concentration and monomer valency 
of the ith layer, respectively. The two above equations 
are solved iteratively for the ith layer concentration 0^, 
while assuming that monomer concentrations from all pre- 
viously adsorbed layers, i — l,i — 2, . . . ,1 are fixed and 
known from previous iterations. These SR interactions are 
contained in the two sums appearing in the right hand side 
of Eqs. lO and Q, Ss and So- The sum 

j=iA-2... 

is the monomer concentration at the point x, summed over 
all similarly charged layers: j = — 2,..., which repel 
the monomers of the ith layer. Similarly, the other sum: 

J2 ^'i^) (6) 

J— i — 3... 

is summed over all oppositely charged layers having an 
attractive SR interaction with the newly adsorbing PEs. 
It should be noted that Ss and So are both functions of 
the layer number i, but the subscript i is omitted for sim- 
plicity. In the following, we consider only monovalent PEs 
and set the odd layers as negatively charged, Z2i+i = — 1, 
whereas the even ones as positively charged, Z2i — 1. We 
also note that for high monomer concentrations the terms 
5*0 and Ss are large, and thus higher orders of both the 
excluded volume and the attractive interactions may be 
necessary. However, in this simple model we restrict our- 
selves to the third order only. 

The solution of the pair of 2nd order differential equa- 
tions, Eqs. ||2J) and requires four boundary condi- 
tions. Two of them are for the bulk where we choose 
4>i{x — > 00) = 0, corresponding to a negligible amount of 
PE in the bulk solution (dilute solution), and zero value 
for the electrostatic potential, ({x 00) = 0. At the solid 
surface, a; = 0, we use the electrostatic boundary condition 
dC/d^l^^Q = —AttIbct, where a is the surface charge den- 
sity. For the first PE layer that adsorbs directly onto the 
solid surface {i = 1), we impose the Cahn-de Gennes at- 
tractive boundary condition piS^ dln((/))/dx|^^Q — —d~^, 
where d is a characteristic length for the SR interactions 



between the surface and the PE chains, and d > corre- 
sponds to an attractive surface. 

For all subsequent layers, i > 2, we expect the PE 
chains to partially penetrate into the previous layers be- 
cause of the complexation. In order to avoid the possibility 
of fully interpenetrated layers, we introduce a hard wall 
for each layer at an arbitrary location x*. Otherwise, in 
our case the cationic and anionic polyelectrolytes would 
form a neutral complex, which will not be able to form 
a stable PE multilayer. Because it is known from experi- 
ments that the adsorbing PE chains of any specific layer 
do not fully mix with the previous layers we introduce the 
concept of the hard wall. Its justification would require 
further studies. For the ith deposited layer, there is a (ar- 
tificial) hard wall x* inside the previous layers so that no 
monomers from the ith layer can reach the region x < x* . 
In order to simplify notation, the layer index i is dropped 
from the hard wall notation, x* . 

The adsorption of every layer brings about a reversing 
of the overall charge of the surface-PE-small ion complex. 
When the adsorbed amount of ions and PE chains exactly 
balances the total charge of the surface and previous ad- 
sorbed layers, the electrostatic field perpendicular to the 
surface is exactly zero (Gauss law). The hard wall of the 
adsorbing PE layer x* is taken somewhat arbitrarily as 
the point where the electric field is zero. As an example, 
the location of x* is depicted in Fig. 2. 

Our choice of x* is motivated by our understanding 
the complexation procedure. The driving forces for the ad- 
sorbing PEs to penetrate the preceding layer are the elec- 
trostatic attraction between oppositely charged PE chains 
and the ion pairing between charged monomers. This elec- 
trostatic attraction is driven by an attractive electric field 
for X > X* . For x < x* , the electrostatic field repels the 
adsorbing PE, and no significant complexation in that re- 
gion is expected. Since we assume that the SR attraction 
is a result of non-equilibrium complexation between elec- 
trostatically attracted PE chains, we do not take x* to be 
dependent on x- It is important to note that such complex- 
ation cannot occur between similarly charged PE chains, 
because the repulsive interaction between the two chains 
as well as the excluded volume repulsion would drive the 
chains to separate rather than inter-penetrate. 

The numerical procedure used to solve Eqs. Q said 
as applied to a single adsorbing PE layer, is based on 
the relaxation method (54j , and was presented in detail in 
a previous publication ^Hl- Here, we use the same proce- 
dure for the layer-by-layer build-up. After obtaining the 
solution for the first PE layer concentration, 4>f{x), and its 
resulting potential, Cix)^ the layer monomer concentration 
profile 0^ (x) is frozen and added as a charge density source 
to the right-hand-side of Eqs. (|3I4|) . These equations are 
now solved for the second layer using the hard wall x* 
and the two virial coefficients, v and x- The procedure is 
then repeated iteratively for all following layers in order 
to obtain a multilayer stack. 
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Fig. 2. The electrostatic potential ( (in dimensionless units) 
of the first four alternating PE layers is presented as a function 
of the distance from the surface x. In each extremum point, the 
electrostatic field d(^/dx changes sign, allowing the multilayer 
to attract an oppositely charged polyelectrolyte. The location 
of the farthest peak (marked by a dashed line) is taken as 
x'^ = X* — the hard wall condition for the interpenetration 
of the next (fifth) layer. This electrostatic potential was taken 
from the numerical profiles such as in Fig. 3. All parameter 
values are specified in Fig. 3. We note that the layers shown 
here are the initially adsorbed layers, and hence the strength of 
the potential oscillations is strongly affected by the existence 
of a charged wall. The amplitude of the electrostatic potential 
oscillations in the more distal region increases towards the final 
layers, due to the increase in the adsorbed layer amount and 
charge. 



3 Results 



CO 



CO 




800 



Fig. 3. The formation of a multilayer stack is shown for high x 
value, X = 0.55a^, and the choice of parameters: Cgait = l.OM, 
a = lOA, o- = 2 • 10"^ A"^, d = 50A. For clarity purpose, 
we show only three groups of layers: layers 1-16 (left), 29-33 
(middle) and 46-50 (right). For both polymers v = O.OSa'^, 
uj'^ = 0.5a®, / — 0.5. The aqueous solution has e = 80 and 
T — 300K. The polycation profiles are marked by a solid 
line, while those of the polyanions by a dashed line. Three 
regions can be seen in the graph. Near the surface the mul- 
tilayer concentration decays rapidly for layers 1-5 (proximity 
region), and then increases rapidly in layers 6-10 (intermediate 
region). The third region is where the multilayer concentration 
stabilizes. This stabilization occurs at a higher value than in 
the initially adsorbed layers. The layers in the distal region 
are highly interpenetrating, so that any layer interacts with 
about five other layers during its adsorption process. Note that 
the lowest monomer volume fraction in the proximity region is 
0.025, which is much lower than in the distal region. However, 
this layer is strongly complexated with the previous layers and 
should still be dense enough to survive the washing procedure. 



Our calculations show a strong dependence of the mul- 
tilayer formation on the value of the SR attraction coef- 
ficient X ill the case of a weakly good solvent, modeled 
via the 2"^^ virial coefficient v — O.OSa"^ where a is the 
monomer size. For low amounts of added salt Cgait = O.IM, 
the formation of multilayers requires very large x/'^^ ^ 3 
values, while for higher amounts of salt Csait = IM the 
required x values drop to more realistic values of x/o^ ~ 
0.4 — 1. For all salt concentrations, low % values cause 
the adsorbed amount of monomers in each layer to decay 
strongly with the layer number, so that very few layers 
are formed. For high x values, the amount of adsorbed 
monomers decreases for the initial layers and then in- 
creases back. The adsorbed amount in each layer is found 
to reach a stable value because of the third virial term. 
The threshold x value is shown below to depend on the 
amount of salt in the solution as well as the initial surface 
charge and the monomer charged fraction. 

The numerical solution of Eqs. JSJ and l@J yields the 
formation of multilayers, as presented in Fig. 3, under the 
proper choice of parameters. As can be seen from the fig- 



ure, the multilayer can be divided into three spatial re- 
gions. In the proximity region, containing the first few 
layers, the adsorbed amount decreases substantially. In 
the intermediate region (layers 6-10), the monomer con- 
centration increases rapidly to much higher values. In the 
distal region, (under some conditions discussed below) the 
adsorbed amount stabilizes, and the multilayer formation 
continues. The adsorbed layers are shown to be very wide 
(of the order of tens of nanometers) and highly concen- 
trated. The interpenetration between the layers looks to 
be quite significant. The location x* , where the next layer 
begins to adsorb, is shared by monomers from all four 
previous layers. This strong interpenetration is the driv- 
ing force of the multilayer formation, since it allows for a 
strong interlayer SR attraction. Without it no significant 
overcharging is achieved. The overall charge of each ad- 
sorbed layer is much higher than the initial surface charge, 
showing that there is no exact charge reversal in PE ad- 
sorption. 
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Fig. 4. The multilayer profile with lower x value, x = 0.39a , 
and all other parameters as in Fig. 3. For this low x, the 
decrease in the non-electrostatic interaction between the PE 
chains causes the adsorbed layers to decay rapidly after four 
layers, and no stable multilayer stack is formed. 



The multilayer formation is characteristic of high x 
values. In the opposite limit of low % values, no stable 
multilayer stack is formed because the complexation be- 
tween the layers is not strong enough. This case is shown 
in Fig. 4, which is obtained for similar parameters as Fig. 3 
except for a lower SR interaction coefficient The figure 
shows that the adsorbed amount in each subsequential 
layer decays rapidly, until an additional layer cannot be 
adsorbed, and the formation of a stable multilayer stack 
is not possible. 

Within our model the formation of a stable multilay- 
ered stack requires a third order virial term. Only when a 
strong enough third-virial coefficient, of order uj^ ^ 0.5a^, 
is added to the SR interaction term, the multilayer concen- 
tration stabilizes at high, but still physical, values. When 
the third virial coefficient is too low, the adsorbed layer 
concentration does not saturate, and rather reaches un- 
realistic high values. Our calculations also show that an 
increase in the second virial coefficient is not enough to 
stabilize the multilayers. It just drives up the threshold 
value of X- The spatial region where the adsorbed amount 
stabilizes is the multilayer distal region, and it can be 
continued for as many as 80 layers (in our calculations) 
without any noticeable decay in the adsorbed amount in 
each layer. 

We end this section by showing three further results. 
In Fig. 5 we show the overall thickness of the adsorbed 
layers from Fig. 3 as a function of layer number. In the 
mean-field model, the thickness of the adsorbed layer is 
taken as the position of the last monomer concentration 
peak, which is a lower estimate for the layer width. The ad- 
sorbed layer width increases weakly for the first few layers 
(proximity and intermediate regions), and then increases 
almost linearly with the layer number (distal region) for 
the entire 50 layers. The linear increase shows that the 
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Fig. 5. The total width D of the adsorbed multilayer in Fig. 3 
is plotted as a function of the layer number. The (incremental) 
layer width is extracted from the peak position in the monomer 
concentration. The total width D is seen to increase mildly 
for the initially adsorbed layers (proximity and intermediate 
regions), and then increase almost linearly, corresponding to 
the stable multilayer formation, with a constant thickness per 
each adsorbed layer. 



multilayers are indeed stable and reaches very high layer 
numbers. 

In Fig. 6 we present the threshold strength of x that 
is needed to form multilayers as a function of the added 
salt amount. As can be seen from this graph, an increase 
in the amount of added salt causes the necessary x to de- 
crease strongly for low salt concentrations. For higher salt 
concentrations the x value is almost constant. The depen- 



dence of X on Csait, X 



^'salt 



fits roughly a power law with 



a ~ —0.8, as can be seen in the inset of Fig. 6. However, 
this empirical scaling is valid only for small range of salt 
concentrations. 

The dependence of the threshold x value on the mono- 
mer charged fraction / is presented in Fig. 7. The x thresh- 
old increases with the increase of /. Here, too, the nu- 
merical results do not imply any simple scaling relation 
between / and x (see inset of Fig. 7). 



4 Discussion 

The spatial behavior in Figs. 3 and 4 can be explained 
by the following argument. In the proximity region, close 
to the surface, the adsorbed layers have high monomer 
concentration and small width. This small width does not 
allow for significant complexation between the adjacent 
layers. Therefore, it causes the SR attraction between the 
already adsorbed layers and the adsorbing PE chains of 
the current layer to be low. This, in turn, causes a decrease 
in the concentration of adsorbing monomers, accompanied 
by an increase in the layer width. 
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Fig. 6. The threshold value of x needed to create stable multi- 
layers is plotted as a function of Caait • For low salt we get high 
values oi X ^ Sa"^, while for high salt the value of x has lower 
values, around 0.5a'^. All other parameters are as in Fig. 3, 
except d = lOA. The inset shows the same dependence on a 
log-log plot. For low Csait values, the slope of the line can be fit 
to X ~ , but for higher salt concentrations the exponent 
becomes lower. Since these changes occur over a single decade 
in Csait values, there does not appear to be a good scaling law 
for X function of Csait. 
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Fig. 7. The threshold value of x needed to create stable mul- 
tilayers is plotted as a function of /. For low / values, the 
monomer-monomer repulsion is low, and only a weak short- 
range attraction between the polymers is needed. For higher / 
values the threshold value of x increases. All other parameters 
are as in Fig. 3, except d = lOA. The inset shows the same de- 
pendence on a log- log plot. No clear scaling law can be found 
for the dependence of the threshold of x on /. 



The behavior of farther layers depends crucially on 
the strength of the short-range attraction, x- For low x 
values, the monomer concentration in the farther layers 
continues to decay, and the conditions are insufficient for 
multilayer formation. This situation is shown in Fig. 4. 
In the opposite case of large enough x values (depicted 
in Fig. 3), the increase of the layer width causes the ad- 
sorbing polymers to interact with more than one adsorbed 
layer. The complexation between adjacent layers becomes 
stronger, allowing the monomer concentration in the ad- 
sorbing layer to increase beyond that of previous layers 
(see intermediate region of growth, layers 6-10 in Fig. 3). 
When the monomer concentration in the adsorbing layer 
becomes high enough (layers 10 and above in Fig. 3), the 
SR attraction between the different PE chains is balanced 
by the third virial term of the excluded volume, and the 
adsorbed amount in each additional layer stabilizes. This 
stable multilayer is characteristic of the distal region, and 
persists to dozens and even hundreds of layers without 
any noticeable decay. We note that the specific built-up 
of the first dozen layers is a direct consequence of our 
simple model of attractive interactions. A more elaborate 
model may be needed for quantitative comparison with 
experimental findings. 

We now turn to the x threshold value needed to obtain 
a stable multilayer formation. The threshold comes about 
because of the competition between the SR and electro- 
static interactions. The decrease of the x threshold with 
salt, as seen in Fig. 6, can be explained in the following 
way. Added salt screens the electrostatic interactions and 
results in an increase in the PE adsorbed amount and layer 
width The thicker layers have a larger contribution 
to the attractive SR term in Eq. ^ and lead to a lower x 
threshold. 

The dependence of the threshold % on / (Fig. 7) is 
explained qualitatively as follows. An increase in / causes 
an increase in the monomer-monomer repulsion between 
the adsorbing PE chains, and an increase in the electro- 
static attraction of the adsorbing PE chains to the al- 
ready adsorbed oppositely charged layer. The increase in 
the threshold x with / shows that the main effect of in- 
creasing / is to decrease the adsorption, meaning that 
the main driving force of multilayer formation is not the 
charge reversal caused by the adsorbing polymer layers, 
but rather the SR interaction. It is important to note that 
in experiment 40 a threshold / value for the multilayer 
formation was found. Above this threshold, the multilayer 
concentration decreases with /, which is in agreement with 
our findings. A threshold in the / value is not found in our 
calculations, mainly because our SR interactions are ex- 
ternally imposed and do not depend on the value of /. 
However, the threshold can be understood qualitatively 
as the value of / for which the polymer chains begin to 
interpenetrate, giving rise to the SR attraction between 
them. We believe that further studies are needed to bet- 
ter understand the origin of this threshold value. 

Our simple model is subject to several limitations. 
First, we use the mean- field theory and the ground-state 
dominance approximation, valid for long PE chains. The 
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adsorption of short polyelectrolyte chains requires other 
treatments such as molecular dynamics, Monte Carlo sim- 
ulations or lattice models. Second, during the adsorption 
of each layer, we assume that all preceding layers are 
frozen, meaning that they do not dissolve back into the so- 
lution or change their spatial conformation. This assump- 
tion stems from the fact that within mean-field theory 
there is no way to distinguish between an adsorbed chain 
and a chain that is merely 'stuck' at the surface vicin- 
ity. This deficiency of mean-field does not allow us to give 
an accurate model for the washing procedure. Experimen- 
tally, we expect the polymer concentration to decrease, es- 
pecially during the washing step, as was modeled by other 
techniques However, due to long relaxation times we 
do not think that the washing will affect drastically the 
structure of the already adsorbed multilayer stack. Finally, 
we use a very simple model for the SR attraction, and do 
not offer any explanation for its dependence on the PE 
parameters. Despite these limitations, we believe that our 
model gives an insight for the multilayer formation prob- 
lem. 

Quite recently, we have become aware of a alterna- 
tive mechanism for multilayer formation suggested by Q. 
Wang (55) . This involves multilayer formation for poor sol- 
vent condition (balanced by electrostatics). The multilayer 
formation can be achieved even if the polyanion and poly- 
cation chains repel each other at short distances, because 
the poor solvent condition induces stable multilayer built- 
up. The solvent condition in Wang's model plays a similar 
role as the attractive interaction between the cationic and 
anionic chains in our model. 



5 Conclusions 

We present a model aiming to explain PE multilayer for- 
mation for marginally good solvents. The model is based 
on strong enough short-range interactions between the 
polyanion and polycation. This strong short-range interac- 
tion is shown to be indispensable for our modeling of such 
stable multilayers. We show that the multilayers form eas- 
ily in high ionic strength conditions, and that their forma- 
tion does not rely exclusively on the electrostatic attrac- 
tion to the previously adsorbed layers. We also calculate 
what is the threshold strength of the short-range interac- 
tions needed for the formation of multilayers as a function 
of the salinity as well as monomer charge fraction. 

In our model the multilayers are quite thick and in- 
terpenetrating, while in the experiments the layers are 
thinner. However, we believe that this simple model gives 
good insight on the problem of multilayer formation, and 
can serve as a starting point for more refined models. A 
possible extension will be to use a more specific model for 
the short-range interactions between the PE chains, which 
may give a better explanation to the experimentally ob- 
served multilayer formation. 
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